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Abstract 

Solid state qubits promise the great advantage of being naturally scalable to large quantum 
computer architectures, but they also possess the significant disadvantage of being intrinsically 
exposed to many sources of noise in the macroscopic solid-state environment. With suitably chosen 
systems such as superconductors, many of sources of noise can be suppressed. However, imprecision 
in nanofabrication will inevitably induce defects and disorder, such as charged impurities in the 
device material or substrate. Such defects generically produce telegraph noise and can hence be 
modelled as bistable fluctuators. We demonstrate the possibility of the active suppression of such 
telegraph noise by bang-bang control through an exhaustive study of a qubit coupled to a single 
bistable fluctuator. We use a stochastic Schrodinger equation, which is solved both numerically 
and analytically. The resulting dynamics can be visualized as diffusion of a spin vector on the 
Bloch sphere. 

We find that bang-bang control suppresses the effect of a bistable fluctuator by a factor roughly 
equalling the ratio of the bang-bang period and the typical fluctuator period. Therefore, we show 
the bang-bang protocol works essentially as a high pass filter on the spectrum of such telegraph 
noise sources. This suggests how the influence of l//-noise ubiquitous to the solid state world 
could be reduced, as it is typically generated by an ensemble of bistable fluctuators. Finally, we 
develop random walk models that estimate the level of noise suppression resulting from imperfect 
bang-bang operations, such as those that cannot be treated as (5-function impulses and those that 
have phase and axis errors. 

PACS numbers: 03.65.Yz, 03.67.Lx, 05.40.-a 
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I. INTRODUCTION 



In order to implement solid-state quantum information processing devices, the noise 
sources causing decoherence of their quantum states have to be carefully understood, con- 
trolled, and eliminated. This is a formidable task, as a solid-state environment generically 
couples a macroscopic number of degrees of freedom to any such device. Thus, a funda- 
mental prerequisite for any design is that it must significantly decouple the quantum states 
used for computation from phonons and other quasiparticles in the underlying solid crys- 
tal. Examples of such designs are those employing discrete states in quantum dots [l| or 
superconductors with a gapped density of states P, S, Q 0. Q 

Most research going beyond this fundamental prerequisite has concentrated on decoupling 
devices from external noise sources such as electromagnetic noise generated by control and 
measurement apparatus [7|. On the other hand, there inevitably are internal noise sources 
because the fabrication of gates, tunnel junctions, and other functional components cre- 
ates defects in the underlying crystal. Prominent examples of such defects are background 



charges in charge-based devices or cricital current fluctuations in flux-based devices 
A clear signature of such defects is telegraph noise in the case of a few defects or 1/ /-noise 
in the case of a larger ensemble ^(|. With the growing success in engineering the electro- 
magnetic environment, these defects are becoming more and more the key limiting sources 
of decoherence. 

Such defects do not fall in the large class of noise sources that can be approximated well as 
a bosonic bath, and this fact complicates analysis. Whereas it is realistic to treat a bosonic 
bath in the tractable near-equilibrium thermodynamic limit where fluctuations are purely 
Gaussian [ill . 12 . kJ, localized noise sources with bounded spectra like the defects in which 
we are interested produce noise that is significantly non-Gaussian. Theories treating large 



I l|, Q|. However, with the 



ensembles of non-Gaussian noise sources have been presented 
ongoing improvement in nanofabrication technology, it is realistic to consider the case where 
non-Gaussian noise sources are reduced down to only a single one or a few per device. This 
is the case we treat here, and thus the defects find a more realistic representation as a small 
set of bistable fluctuators [16] (henceforth abbreviated bfls). In principle, this approach 
can be extended to larger sets of bfls with a range of different mean switching times (e.g., 
an ensemble with an exponential distribution of switching times that produces l//-noise 
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3). 

This report is organized as follows. Section II presents the model of a single bfl in 
the semiclassical limit, where it acts as a source of telegraph noise. Section III introduces 



an idealized open loop quantum control technique, quantum bang-bang control 
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221 ] . which is suitable for slowly fluctuating noise sources. Section IV explains how we 
simulated the qubit dynamics under the influence of noise with and without bang-bang 
control by integrating of the corresponding time-dependent Schrodinger equation. As a 
measure of the decoherence, we analyze the deviations of the qubit's trajectory on the Bloch 
sphere from that of the noiseless case. These deviations take the form of a random walk 
around the noiseless-case trajectory. We therefore analyze the suppression of these deviations 
by comparing the variances of these random walks with and without bang-bang control. 
Both numerical and analytical solutions (the latter in the long-time or "diffusion" limit) 
are presented. Comparison of the numerical simulations to the analytical solutions shows 
excellent agreement. We then analyze how these results change when practical limitations 
are considered such as the fact that a bang-bang pulse cannot be an ideal 5-function impulse 
and the fact that the duration or polarization axis of the pulse may suffer from random 
fluctuations. We show at the end of Section IV. B that within large margins bang-bang 
suppression of the bfl noise is not inhibited by having a finite, rather than infinitesimal, 
pulse length. However, in Section V, we do find that duration and polarization axis errors in 
the bang-bang pulses can significantly affect the suppression of bfl noise. We present a point 
of optimum performance. Section VI concludes with remarks on several recent publications 
concerning the suppression of telegraph or 1/ /-noise. 

II. MODEL OF THE BISTABLE FLUCTUATOR IN ITS SEMICLASSICAL LIMIT 

We describe the bfl-noise influenced evolution of the qubit in its semiclassical limit by 



using a stochastic Schrodinger equation 
nian 

reff 



23, 



with the time-dependent effective Hamilto- 



Hf(t) = H q + H noise (t) (1) 
H q = he q a? + HA q a« (2) 

#noise(t) = ha £bfl(0 (3) 
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where e q and A q define the free (noiseless) qubit dynamics. £bfl(£) denotes a function 
randomly switching between ±1 (see Fig. |TJ), which represents a telegraph noise signal. 
The switching events follow a symmetrical Poisson process, i.e., the probabilities of the 
bfl switching from +1 to —1 or —1 to +1 are the same and equal in time. The Poisson 
process is characterized by the mean time separation t^q between two bfl flips. The coupling 
amplitude to the qubit in frequency units is a. The relation of this Hamiltonian to a 
microscopic model is explained in the Appendix. 
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FIG. 1: Schematic plot of a typical Poisonian bfl noise signal and its resulting random walk behavior 
(in the limit of small deviations). The periodic fast switching step function represents a bang-bang 
pulse with a time scale ratio: t^/t^ = 10 and yields a quite smaller random walk step-length. 



r Sys 



denotes the evolution period of the qubit in the noiseless case. 



Starting with an arbitrary initial state of the qubit, represented by some given point 
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on the Bloch sphere, we can numerically integrate the corresponding stochastic differential 
equation and obtain the corresponding random walk on the Bloch sphere 

&{t) = Texp^-i/fr^ Hf(s)ds S ja(0) (4) 

with T denoting the usual time-ordering operator. 




FIG. 2: Schematic plot of noisy qubit evolution generated by Poissonian telegraph noise. The 
resulting random walk (dot-dashed line) on the Bloch sphere is comprised both of deviations 
Aod^ in parallel to the free precession trajectory (dotted line), which correspond to dephasing, 
and deviations Acr re i perpendicular to it, which correspond to relaxation/excitation. 
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III. BANG-BANG CONTROL PROTOCOL 



We propose to reduce the influence of the bfl-noise by applying to the qubit a continuous 
train of 7r-pulses along the <r x -axis. This refocusing pulse scheme essentially corresponds to 



the standard quantum bang-bang procedure 20, 2.1, 22] or the Carr-Purcell-Gill-Meiboom 



echo technique from NMR 25|. For technical convenience, we consider the 7r-pulses to 
be of infinitesimal duration. This simplification is not crucial as will be detailed later in 
Section IV.B. The pulses are assumed to be separated by a constant time interval Tbb- The 
mean separation t\,q between two bfl flips is assumed to be much longer than Tbb- For 
theoretical convenience, we also assume that im is shorter than the free precession period 
of the qubit. This too is not a crucial restriction. (It can always be overcome by changing 
to a co-precessing frame.) 

Qualitatively, bang-bang control works as follows. Since r b b <C r h{i , it is usually the 
case that the bfl does not flip during the time between two bang-bang pulses that flip the 
qubit. In this way, the bang-bang pulses average out the influence of H noise (t). In fact, 
the refocusing scheme fully suppresses the cr 2 -term of the static Hamiltonian (J21 (compare 
Fig. EJ); but this turns out to be no crucial obstacle to universal quantum computation as 
will be outlined later on. As one can visualize in Fig. ^ it is only when a bfl flip occurs 
during a bang-bang period that the net influence of the bfl felt by the qubit is nonzero, 
and the qubit thus suffers some random deviation from its trajectory in the noiseless case. 
Taken together, these random deviations constitute a random walk around the noiseless case 
trajectory. While this walk is actually continuous, it can be modelled as a discrete walk with 
steps that are randomly distributed in time, one step for each bfl flip (see e.g. The 
average step length is essentially the product of the noise coupling strength a and the mean 
time the bfl in its present state can influence the qubit. Without bang-bang control, this 
mean influence time is Tbfl, whereas with bang-bang control, it is reduced to Tbb- Therefore, 
both with and without bang-bang control, the random walk has the same time distribution 
of steps, but with bang-bang control the step size can be significantly reduced roughly by a 
factor of the ratio of time scales r hh /r h{i . 
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IV. RANDOM WALK ON THE BLOCH SPHERE 



Now we study this proposal quantitatively. We simulate these random walks both 
with and without bang-bang control by integrating both numerically and analytically the 
Schrodinger equation, Equ. (4), with the stochastic Hamiltonian of Equs. (1-3). As generic 
conditions for the qubit dynamics, we choose e q = A q = Q . Without loss of generality, we 
set the qubit's initial state to be spin- up along the z-axis. If the qubit-bfl coupling a were 
zero, then the qubit would simply precess freely on the Bloch sphere around the rotation 
axis a% + &2 (the dotted line in Fig. |2J. Hence, we expect for a sufficiently small coupling 
(a <C Qq) only a slight deviation of the individual time evolution compared to the free evo- 
lution case (the dashed line in Fig. |2J). For the coupling strength, we take a = 0.1f2 . All 
the following times and energies are given in units of the unperturbed system Hamiltonian, 
i.e., our time unit r Sys is given according the free precession time 7iT Sys /\/2, and our energy 
unit is given by AE = y / e q H-A^ = y/2D, . The time scale ratio is taken to be r bfl /r bb = 10 
if not denoted otherwise. 

This approach accounts for the essential features of our specific situation: the long correla- 
tion time of the external noise, essentially T b fl, its non-Gaussian statistics and its potentially 
large amplitude at low frequencies. These properties are crucial and are difficult, although 
not impossible, to take into account in standard master equation methods. 

A. Numerical simulations 

We have numerically integrated Equ. (J3J) and averaged the deviations of the random 
walk evolution from the unperturbed trajectories for times up to 100rs ys over N = 10 3 
realizations. Larger simulations have proven that convergence is already sufficient at this 
stage. We shall examine the root-mean-square (rms) deviations of this ensemble at given 
time points 



with and without bang-bang control. In other approaches, such as those based on master 
equations, one separates dephasing and relaxation. Both are contained here in Equ. (0). We 
shall point out notable differences between these two channels. The deviation as a function 
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FIG. 3: Time evolution of the rms deviations for bfi- induced random walks with and without 
bang-bang control at a coupling constant a = 0.1 and a typical flipping time scale Tbfl = 10 _2 rs ys . 
The separation between two bang-bang pulses is Tbb = 10 _3 rs ys . The straight lines are square- 
root fits of the analytical derived random walk model variances (plotted as triangles). Inset: 
Components of the deviations from the free precession trajectory that are parallel to it (dephasing) 
and perpendicular to it (relaxation/excitation) with bang-bang control. 



The total deviations on intermediate time scales are suppressed by a ratio on the order of 
10. A detailed numerical analysis shows that without bang-bang suppression, the deviations 
parallel to the free precession trajectory (which correspond to dephasing) are of similar size 
to those perpendicular to free precession (which correspond to relaxation/excitation). In 
contrast, with bang-bang control, dephasing is almost totally absent as one can see in the 
inset of Fig. El 

The main double-logarithmic plot of Fig. 3 shows that on short time scales (t < 0.1r Sys , 



which corresponds to < 10 random walk steps), deviations increase almost linearly in time. 
It is not until times on the order of rs ys that the noise-induced deviations start to behave as 
typical classical random walks, increasing as a square-root in time. 



B. Analytical random walk models 

We now develop analytical random walk models for our system. The random walk on the 
Bloch sphere is in general two-dimensional, consisting of both parallel and perpendicular 
deviations to the free evolution trajectory. Bang-bang control, as was seen in the above 
numerical results and as will be seen in the following analytical results, essentially reduces 
the random walk to one-dimension as only the perpendicular deviations remain significant. 
In the following, we restrict ourselves to the long-time (many random walk steps) regime. 

We first calculate for both cases the probability distributions of the deviations after one 
bfl flip ("one-step deviations" in terms of the discrete random walk). The fluctuation of the 
period between r^, r leads to dephasing, which can be evaluated at a <C e q , A q to 

A^deph = 2tt cos -) rm * ±2-^L_ aTbfl ( 6 ) 

V r pcr T per / ^ q ~r £ q 

where the prefactor cos0 = , A „ q = takes the effective trajectory radius into account. 

For the relaxation/excitation effect of the noise, one has to use the projection of the 
perturbation orthogonal to the free axis, using sin?] = 2 °^ 2 . Furthermore this type of 
deviation also depends on the actual position of the spin on the Bloch sphere, e.g. there is 
no relaxation when the state is at one of the poles. 

Averaging in rms-fashion over a full azimuthal cycle leads to a factor of l/y/2. Moreover, 
the impact of relaxation/excitation is scaled down by an additional factor of cos0 = , Aq 

V A q+ e q 

corresponding to the projection of the Bloch vector onto the precession axis, which fur- 

'per — 'per 



thermore decreases the deviation angle. In total, using ri. ~ r per to first order in a, we 



find 

1 T fl A 3 

Aa™ = 2tt cos sin r]— = cos <p^§- ~ V2 2 ? 2 \3/2 aThR - ( 7 ) 

V 2 Tper l e q + ^qJ 

Our root-mean-square measure of the impact of the noise, Equ. (5), does not handle these 
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FIG. 4: Plot of a typical one-step deviation from the unperturbed qubit trajectory with generic 
values for e q and A q . The fractions of the bfl fluctuations in <r z -direction have to be distinguished 
with respect to their effects on the qubit: those that yield dephasing deviations that are parallel to 
the free precession trajectory (proportional to sinc/>) versus relaxation/excitation deviations that 
are perpendicular (proportional to sinry). Both parts are additionally domineered by a factor of 
cos^ due to the diminished radius of the trajectory starting from the initial state a z = +1. The 
impact of the relaxation/excitation generating part is furthermore depending on cos^ as well as 
sinx, the azimuth angle of the qubits present position. 

two kinds of deviations separately, but rather adds them up to: 



total 



A -bfl 2 , A =» 



bfl 2 



A 2 e 2 
q q 



+ 2- 



A 6 

q 



ar bfl 



1)3/2 



4(A 2 + e 2 ) A 2 e 2 + 2A6«r bfl (8) 
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Our rms treatment disregard the different types of decoherence, dephasing and relax- 
ation/excitation, corresponding to phase and bit-flip errors respectively. This is no crucial 
drawback but merely lies in the nature of our generic situation. If needed, both components 
can be isolated. 

The derivation of the maximal one-step deviation for the bang-bang controlled situation 
has to be handled differently. The deviation resulting from a bfl flip during a bang-bang 
pulse period is maximal if the step happens exactly at the moment of the second qubit 
spin-flip (i.e., in the middle of the bang-bang cycle). When this happens, the refocusing 
evolution has in its first half a drift, for example, to the "right" (compare to Fig. EJ) and in 
the last half an equal aberration. 
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FIG. 5: Sketch of a maximal one-step deviation during a bang-bang modulated cycle, which 
appears if the bfl state flips precisely at the intermediate bang-bang pulse time. The dephasing 
part of deviation evidently averages out, while a relaxating aberrance arise proportional to the 
noise-coupling constant a. 
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The resulting one-step deviation appears to be on the order of 2art,b- However, this is 
scaled down by a factor of 1/ y/2, as the impact of the aberration in x-direction is proportional 
to a factor of sinx, where \ denotes the longitudinal angle of the present spin position on 
the Bloch sphere (see Fig. EJ). This is because the a z -component of the noisy evolution 
does not influence the qubit, if it is near the a z = ±l-state and its influence is suppressed 
correspondingly in between. As we are mainly interested in mean aberrations after many 
random walk steps, we simply average the maximal one-step deviation over one precession 
period in the usual rms manner to obtain 



= V2ar hh . (9) 

Obviously, this variance only contributes to relaxation. 

In the long time limit, we replace the fluctuating number of random walk steps for a given 
time At of noisy evolution by its expectation value iVbfl = At/r^. This allows us to use the 
number of random walk steps as time parameter j^J. This simplification does not introduce 
significant error, as the relative number variation for At scales as ^p—^ = a/ t^/ At — > in 
our preferred long-time limit. We encounter two different one-step-distributions, depending 
on whether the numeration of the step is an odd or even (corresponding to an "up" or 
"down" state of the bfl). For defmiteness, we assume the bfl is initially in its "upper" state, 
which is of no influence on the long time limit as the memory to the initial state is already 
erased. The step-size distribution of the bfl model in our small deviation regime is given 
from Poisson statistics 

^dVenW = eTX/ J ±X) (10) 

with (3 = ^aTbfl the typical one-step deviation as calculated in Equ. |H1 9{x) denotes the 
Heaviside step function. We neglect the correlations between transverse and perpendicular 
deviations as we expect them to average out in the long-time limit. 

For the bang-bang suppressed random walk, the flipping positions of the bfl-noise sign in 
the bang-bang time-slots are essentially randomly distributed as long as Tbb <C Tbfl- That 
is why we find a constant step-size distribution between zero and a maximum divergence of 
7 = (see Equ. EJ), namely 

•S*-<«) = e{±x)e ^ x) . (id 
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By means of these one-step probability distributions, we are able to calculate via the 
convolution theorem the distributions for 2iVbfl-step random walks. Specifically, they are 
the inverse Fourier transforms of the iV bfl -fold products of the Fourier transforms of the 
two-step distribution 2(| . For the case without bang-bang control, we find 

dk ^-ikx ( ' ^ 



_ T 2vr/? 2Jv bfl \l-2 cos(A;)e- 1 //3 + e" 2 ^ 
whereas for the case with bang-bang control, it is 



) (12) 



dk ^/ [l-oosqT + l)^] ^ (13) 



_ w 2vT7 27V bfl ^ [l _ cos(Jfe)] 

with & and J£~ _1 denoting the discrete Fourier transformation and its inverse, respectively. 

Already for random walk step- numbers on the order of 10, the resulting distributions 
are almost Gaussian. Their standard deviations give the rms deviations of the random 
walk models plotted in Fig. El As in the numerical simulations at long times, they grow 
as a square-root of the number of steps. As one can recognize, the underlying two-step 
distributions in the fc-space (i.e., the functions in the large brackets of Equs. (J 12)) and (|13j) ) 
are symmetric and differentiable around zero such that the above integrals can be evaluated 
analytically using the saddle point approximation (the small parameter is k, which is justified 
at least qualitatively in our bounded variable integral). We find for their variances in real 
space representation 

\/5 

A(x bfl (iV bfl ) = v/AW? = 0VM^-ar bfi (14) 
for the case without bang-bang control and 



Aa bb (iV bfl ) = = J ^-OiT hh (15) 

for the case with it. In the large- iV b fl limit, this model shows excellent agreement with the 
numerical simulations. 

At first sight, treating bang-bang pulses as 5-f unction impulses appears to be an ex- 
traordinarily strong assumption, especially because in a physical implementation, the large 
bandwidth associated with very short pulses could excite other noise sources. However, this 
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^-function impulse approximation is only for technical simplification. In fact, going to the 
other extreme of a wide, continuous pulse of the form sin(— t) would also refocus our bfl- 
noise over the course of its periods. Comparing the two-step deviation distributions arising 
from (^-function impulses versus continuous sine waves, one obtains for the 5-function case 



and for the continuous sine wave case 



|7 — x\ 

7 



0(7 - x)0(7 + x) 



(16) 



7T 7T X 

— + — COS 27T- 
/7 47 \ 7 



" ) 1G~ 



3 / x 

Sin I Z7T — 

7 



9(j - x)9( 1 + x){17) 



These distributions are depicted in Fig. El One recognizes that in fact the distribution 
arising in the continuous sine wave case is narrower (and therefore indicates more effective 
noise suppression) than the 5-function impulse case, with the drawback of leaving less free 
evolution time for coherent operation. 




FIG. 6: Comparison of two-step distributions for the random walks with bang-bang control when 
the bang-bang pulses are taken to be <5-functions (left) versus a continuous sine wave sin(— t) 
(right). For clarity, the y-axis is rescaled to the maximum values of the distributions, while the 
a>axis is given in units of 7. 
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C. Distributions of the random walks deviation 

Beyond predicting the variances of the random walks, our analysis also allows evaluation 
of their full probability distributions. We compare them to numerics with and without 
bang-bang compensation by use of simulations with 10 4 realizations at an evolution time 

— r Sys- The numerical histograms of the deviations with their respective one- and 
two-dimensional Gaussian fits are shown in Fig. 




I Ao I I Ac I 

FIG. 7: Histograms of the deviation from free evolution both without bang-bang control (left) and 
with bang-bang control (right). Also plotted are fits to the expected two- and one-dimensional ran- 
dom walk statistics respectively associated with the uncontrolled and controlled cases. Numerical 
data were collected over 10 4 realizations at a fixed time to = T Sys defined such that Tbfl = 0.01rs ys 
and thus iVbfl = TSysAtfl = 100 steps. (NB: The rr-axis scale of the right graph depicting the 
bang-bang controlled case is 15 times smaller than that of the left graph depicting the uncontrolled 
case.) 

We observe that not only the distribution obtained with bang-bang control is much nar- 
rower than the distribution obtained without it, but also that its shape is qualitatively 
different. The maximum of the bang-bang controlled distribution is at zero error. In con- 
trast, the uncontrolled distribution has its maximum at a finite error |Acr| max « 0.01, and 
it has zero probability of zero error. This reflects the one-dimensional nature of the bang- 
bang controlled random walk in contrast to the two-dimensional nature of the uncontrolled 
random walk. 
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D. Bang-bang control working as a high-pass filter 



In order to measure the degree of noise suppression due to bang-bang control, we define 
the suppression factor S tQ as follows for a given evolution time t 

S t0 (r ha /r hh ) = (18) 

We now systematically study the dependence of S to on r^a/r^ for a constant mean 
bfl switching rate rbfl = 10~ 2 T sys at a fixed evolution time to = r sys . The numerical data 
in Fig. |H1 show that the suppression efficiency is linear in the bang-bang repetition rate, 
Sr sys = ^ r bfi/' r bb- The numerically derived value of the coefficient, /i nume ricai ~ 1-679, is in 
excellent agreement with the analytical result /i an aiyticai = a/5/2 ~ 1.581 from our saddle 
point approximation, Equs. (fT^j) and (fTo"]) . 



This small discrepancy between the numerical and analytical results is due to the fact 
that the analytical calculations neglect correlations between the parallel and perpendicular 
components of the random walk. This leads to an underestimate of the rms-deviation 
AffJ 1 ^ in the case without bang-bang control (compare also to Fig. |7|). Therefore, we have 
quantitatively proved our qualitative intuition: bang-bang control affects the bfl noise signal 
like a high-pass filter, an effect that one of the authors has generally predicted for dynamical 
decoupling techniques j^lj]. 



V. LIMITATIONS DUE TO PULSE INACCURACIES 

Thus far, we have tacitly assumed that one could apply perfect, zero-width 7r-pulses 
along exactly the <j x -axis of the Bloch sphere. We now take into account that the control 
pulses themselves typically will have slight fluctuations in their duration or polarization 
that interfere with the desired refocusing. As already shown at the end of Section IV.B, the 
restriction of pulses to infinitesimal duration can be significantly relaxed. We now investigate 
to what extent the restriction to perfect pulses can be relaxed. 
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FIG. 8: The suppression factor 5t (rbfl/rbb) = A^t>w!°l evaluated for to = TSys as a function of the 
ratio of the mean switching time Tbfl and the bang-bang pulse separation Tbb- 

A. Two generic types of bang-bang inaccuracies 

We essentially analyze two generic types of errors that could occur in the control appa- 
ratus when trying to apply 7r-pulses in <j x -direction. One, the duration of each pulse could 
exhibit fluctuations, resulting in fluctuations in the rotation-angle around the desired value 
of 7T. Two, the polarization axis could suffer from directional deviations around the desired 
value of 6" x . Assuming the statistical independence of each pulse error, we expect for both 
types of imperfections a random-walk-like behavior of increasing deviations compared to 
evolutions with perfect pulses. 
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1. One- dimensional pulse error (dephasing) 



We make the quite general assumption that we may model the one-dimensional phase 
fluctuation of the imperfect bang-bang pulses 4>j(x) as a Gaussian distribution of the pulse 
durations and therefore of the rotation angles around their intended value ir. This as- 
sumption should be valid for many physical situations, e.g., if the inaccuracy is due to 
electromagnetic noise in the pulse generator. The Gaussian is parameterized by its standard 
deviation £0 O (see Fig. EJ). Thus, the corresponding pulse angle aberration of the jth step is 
given by 

1 x2 

cj)] d (x) = -= e (19) 



Having assumed a Gaussian distribution, we can exactly evaluate the distributions of the N- 
step deviation A$a? (which are usually given as iV-fold time-convoluted integrals) as follows 
by use of the convolution theorem 



72 



e 2N < (20) 



with 0j d = &[4>j\ denoting the Fourier transform of 0j d and J^~ _1 denoting the inverse 
Fourier transform. 

Therefore the rms displacement in the random walk increases as a square-root in the 
number N of bang-bang pulses: 5<p n = VNd~4>o- Equivalently, the dephasing grows as 
square-root in time 



5<P(t) = Vt/^^o (21) 
on the time scale of our coarse-graining (which is here given as Tbb)- 



2. Two-dimensional pulse error: dephasing and relaxation/ excitation 

A similar argument works when there are also fluctuations around the desired <r x rotation 
axis. Each individual variation of the axis can be split into two components: (1) 50 perp , which 
is perpendicular to the connecting vector between the <r x -axis and the qubit state a(t) on 
the Bloch sphere, and (2) <50 t an, which is transverse to it (see Fig. fTDj> . To first order, the 
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FIG. 9: Sketch of one-dimensional bang-bang aberration. The variations 5(po of the rotation angle 
around the desired value of n leads to slight deviations in parallel to the permitted dynamical 
direction, thus generating dephasing. 

perpendicular part does not disturb the intended spin- flip j33|. However, the transverse part 
does cause a deviation from the ideal spin-flip in a direction toward or away from the previous 
qubit state. (Therefore, it produces relaxation or excitation, as its effect is orthogonal to 
the free <r x -evolution.) Consequently, in a statistical average we only have to consider l/y/2 
of the typical total mean 5<po of the aberration. The effect of a 7r-rotation around an axis 
tilted by an angle <50 tan is a deviation 2<50 t an from the trajectory of the perfect evolutions; 
thus we receive altogether a deviation on the order of a/250o- 
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FIG. 10: Sketch of two-dimensional bang-bang aberration. To first order, variations 5<fip er p of 
the rotation axis perpendicular to the connection vector between a x and the qubit state (here for 
simplicity: <r z = +1) do not influence the intended spin- flip, whereas the variations 5^>tan along this 
line causes deviations on the Bloch sphere perpendicular to the permitted evolution trajectories 
(therefore producing relaxation or excitation). 

Therefore, we obtain analogously to Equ. (|19|) for each single step distribution 

1 x2 ., 

0f (x) = _ _ e (22) 

and analogously to Equ. (j2TJj) for the deviation after iV steps 

1 x2 

$i d = = - e T ^>. (23) 
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Equivalently, in terms of the time t 

5<j>{t) = y/2t/r hh ■ 6<fH>. 



(24) 



B. Numerical and analytical results 

In the same manner as our previous integrations of a stochastic Schrodinger equation, we 
numerically simulate qubit dynamics under inaccurate pulses. In the first instance, we work 
without bfl-noise to verify our analytical random walk model. Later, we add the bfl-noise 
in order to study the competition between the two sources of error. 

1. Random walk due to inaccurate bang-bang pulses only 

We analyze deviations on the Bloch sphere between the noiseless case trajectories that 
occur when the bang-bang pulses are perfect and those when they are not. As per Equ. (JHJ), 
we calculate the rms-deviation over ensembles of iV = 10 3 realizations. As a representative 
time point, we once again choose t = r Sys . This is because, as explained in the discussion 
surrounding Fig. El this time scale should exhibit neither short-time effects nor near-total 
decoherence. From Equs. (J2T|) and (|24j) . it immediately follows that for the mean deviations 
at to if there are phase errors 



As characteristic values for the mean accuracy of single pulses, we choose 5<fto in the 
range of 10~ 6 to 1CT 4 , which should be technologically feasible. As one can see in the double 
logarithmic plots of Fig. ^2 the numerically determined evolutions follow the analytically 
expected square-root type random walk behavior. 




(25) 



and if there are axis errors 




(26) 
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FIG. 11: Plot of the one- respectively two-dimensional imperfectly bang-bang pulsed evolution. 
Dashed lines are square-root fits of the numerical data, while the solid lines denotes the analytical 
calculations. 

2. Random walk due to both inaccurate bb-pulses and bfl-noise 

We now combine our imperfect bang-bang pulse operations with our former bfl-noise 
signal to discuss the applicability of our control scheme when "realistic" pulse generators 
are used. As before, we calculate the rms deviations at to = T Sys by averaging over 10 3 
realizations. The bfl-parameters are those used previously: a coupling strength a = 0.1 and 
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an average switching time t^a = 0.01rs ys . However, with the aim of determining the optimal 
bang-bang protocol in the presence of pulse imperfections, we now consider different pulse 
separation times Tbb/isys ranging from 1CT 5 to 1CT 2 . 

We assume that the errors induced by the bfl and those induced by the pulse generator 
are statistically independent, and thus we sum together both sets of induced deviations in 
the usual rms- fashion. In comparison to the case of ideal bang-bang pulses, Equ. (|T5j). we 
find here the average total deviations induced by both bfl telegraph noise and imperfect 
bang-bang pulses to be: 



in the one-dimensional case where imperfect pulses only impart phase errors (due to impre- 
cise pulse duration), and 



in the two-dimensional case when imperfect pulses impart both phase and relax- 
ation/excitation errors (due to imprecision in the pulses' polarization axes). 

As Fig. El demonstrates, we observe a very good agreement between our numerical and 
analytical results. Such data make it possible to determine an optimal bang-bang separation 
time r£fo. Specifically, this optimum can be derived by calculating the zero value of the 
derivative of Equs. (J27|) and (J28J) with respect to t\^. We therefore conclude that the 
optimal period between bang-bang pulses is: 




(27) 




(28) 




(29) 



for the one- dimensional case and 




(30) 
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FIG. 12: Plot of the Bloch sphere rms deviations received by one-/two-dimensional inaccurately 
pulsed bang-bang compensation of the typical bfl-perturbation. Dashed lines describe the aber- 
rances for pure faulty bang-bang (i.e. without bfi-noise), respectively the exactly compensated 
bfi-case (see Fig. |3J), while solid lines denotes the deviations calculated by random walk analysis. 

for the two-dimensional case. These optimal times respectively correspond to minimized 
variances at to — T Sys of 

A<t = V2 + 1 1/6 ^ (31) 
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for the one- dimensional case of only imprecise pulse durations and 

A< t = V2-V3 + 2 2 /3 ?^01 y fc (32) 

for the two-dimensional case of both imprecise pulse durations and polarization axes. 

VI. CONCLUSION AND REMARKS 

We have investigated the qubit errors that arise from the noise generated by a single 
bistable ffuctuator (bfl) in its semiclassical limit, where it behaves as a telegraph noise source. 
We numerically integrated a corresponding stochastic Schrodinger equation, Equ. (jlj), as 
well as analytically solved (in the long-time limit) appropriate random walk models. As a 
characteristic measure of the resulting dephasing and relaxation effects, we used the rms 
deviation of noisy evolutions compared to noiseless ones. To suppress the effects of this 
noise, we presented a bang-bang pulse sequence analogous to the familiar spin-echo method. 
We claimed this pulse sequence to be capable of refocusing most of the bfl-noise induced 
aberrations. Both in the case without bang-bang control and the case with it, there was ex- 
cellent agreement between our numerical and analytical results on the relevant intermediate 
time scales (i.e., times after a short initial phase where deviations grow linearly instead of 
as a square-root in time, but before the qubit becomes totally decohered). 

In particular, we confirmed our preliminary qualitative picture that bang-bang control 
works as a high pass filter, suppressing qubit errors by a factor /rr bfl /r bb that is directly 
proportional to the ratio of the mean bfl switching time and the period between bang-bang 
pulses. The numerically and analytically calculated constants of proportionality /^numerical ~ 
1.679 and yU ana i ytica i = a/5/2 ~ 1.581 also matched to good accuracy. These results imply that 
the bang-bang procedure is an appropriate remedy against the 1 / /-noise that often is seen in 
solid-state environments. This is because bang-bang control exhibits maximal suppression 
of bfl telegraph noise, and 1/ / noise generally arises from an ensemble of bfls. Finally, one 
has to be aware that also the static <r z -term of the Hamiltonian is averaged out, and this 
generally reduces the degree of control on the qubit. But this is only a technical constraint, 
as one could imagine interchanging two different types of bang-bang pulses (e.g. along the 
x and the y-axis respectively) to admit corresponding quantum-gate operations. 

We previously presented this basic idea in a short paper [27]. The present work extends 
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that short paper by treating the effects of different types on non-ideal bang-bang pulses. 
Moreover, the analytical random walk method is outlined in much more detail, as this 
method should also be applicable to other problems that are difficult to treat in a master 
equation approach. 

Meanwhile, several other extensions of Ref. [27] have been proposed by other research 



groups. Ref. |l7j includes a larger number of fluctuators, described as semiclassical noise 
sources, but restricts itself to a single spin-echo cycle. Ref. analyzes extensively the 
importance of higher, non-Gaussian cumulants and memory effects and arrives at a number 
of analytical results, but it does not treat the option of refocusing. Ref. [2^ treats a full mi- 
croscopic model and compares different variations of the bang-bang pulse sequence. Ref. [3^ 
also treats a full microscopic model with potentially many fluctuators using a Lindblad-type 
approach and covers a wide range of ratios between the fluctuator and bang-bang pulse time 
scales. One of its main conclusions is that a Zeno effect is found in a parameter regime not 
covered by our work. Note that all of these other extensions of our work treat only the case 
of ideal bang-bang pulses. 



Acknowledgments 

We are indebted to G. Falci, I. Goychuk, T.P. Orlando and J. von Delft for helpful 
discussions. Our special thanks goes to A. Kaeck for pointing out the capability of our 
random walk model to describe deviations due to imperfect pulsing. HG and FKW also 
thank T.P. Orlando for his great hospitality during their stay at MIT. WMK gratefully 
acknowledges fellowship support from the Fannie and John Hertz Foundation. This work 
was supported by a DAAD-NSF travel grant and by the National Security Agency (NSA) 
and Advanced Research and Development Activity (ARDO) under Army Research Office 
(ARO) contract P-43385-PH-QC as well as by the DFG through SFB 631. 



APPENDIX A 



We shall now connect the model of a single bfl as a telegraph noise source to a microscopic 



Hamiltonian. We start with the conventional Hamiltonian model of a single bfl, e.g., 



3ll,Q 



14. 



The original qubit is influenced by noise from another qubit, the bfl, which itself is 
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coupled to a thermal environment by a bilinear spin-boson type interaction: 

H = i/qubit + -^qubit.bfl + #bfl + -^bfi,env + H env . (Al) 

where 

^qubit = he q a* + hA q a* (A2) 

^qubit.bfl = haa^a™ (A3) 

H hR = ht m af + hA ha a™ (A4) 

^bfl.env = £Aa z bfl (s f + %) (A5) 

3 

H eav = hj^ujj [a)a 3 + 11/2) (A6) 

j 

The scalar a denotes the coupling strength between the original qubit and the bfl, while the 
scalar A indicates the influence of the environmental heat bath on the bfl. 

It is not obvious how to treat such a combined open quantum system fl Q- The 
common approach of deriving a master equation for the reduced qubit system does not 
work, as it is not clear how to introduce an open quantum system "bfl" as the environment. 



Gassmann, et al. present four alternative approaches jl4J |. Their first approach is to derive 
a standard Markovian master equation for the combined open system "qubit + bfl" and 
trace out the parameters of the bfl afterwards. Their second approach is to consider the 
qubit as influenced by an effective bfl-bath environment by use of an Markovian and secular 
approximation in the limit of small a. Their third approach, which is both the most general 
and the most complicated, is to deduce a master equation by applying a non-Markovian 
weak-coupling perturbation ansatz in second order in a. 

For our investigations, we prefer this last and most general approach: a stochastic treat- 



ment em 



nploying an appropriate randomly changing bfl-noise Hamiltonian term (compare 
also to 1 171). This choice is not only because of practical reasons (to make our numerics 

n □ 

feasible), but also due to empirical considerations (see [8j, |9|], where characteristics of tele- 
graph noise were observed and attributed to bfls). Hence, we restrict our analysis to the 
limit A > a, ie., the limit where the coupling of the bfl to the external environment is much 
larger than interaction between the bfl and the qubit. For convenience, we assume the bfl 
is in its high-temperature limit. (Note that this does not necessarily mean the qubit is also 
in a high-temperature regime for the qubit 's energy scale might be much larger than that 
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of the bfl.) We thus assume the bfl behaves like a classical (i.e., decohered) noise source, 
and we specifically describe the bfl's influence on the qubit with the following stochastic 
Hamiltonian: 

# qub it,b fl ^ H^ a (t) = ha d* (A7) 

In the equation above, £bfl(£) is a random function of time representing the switching of the 
a\ -value between ±1. In our high-temperature limit, we assume ^bfl(^) has symmetrical 
Poissonian statistics (i.e., the probabilities of the bfl switching from +1 to —1 and from — 1 
to +1 are equal and constant over time). Such a symmetrical random process is readily 
described by just one parameter: the typical time separation r bfl between two bfl flips (see 
Fig. [ID. 

The high temperature limit is not a crucial constraint. Treating the strongly thermally 
coupled bfl in an intermediate temperature regime would only result in some asymmetri- 
cally switching £bfl(^)- The typical switching times time for switching the bfl up and 

down respectively satisfy the detailed balance relation ^f- = e~ SEba ^ kbT , where 5-Ebfl de- 

r bfl 

notes the energy separation of the two bfl-states, and T the temperature of the heat bath 
which drives the switching of the bfl. The microscopic structure of the rates depends on 
details of the experiment. Typically, they will be golden rule rates containing the density 
of states of the heat bath and the matrix element of its coupling to the bfl. If that bath is 
made of harmonic oscillators with an ohmic spectral density, we e.g. expect switching rater 
( r bflJ = ^ e ±sE^k^T _^ where «o is the dimensionless coupling strength to the Ohmic 
bath. This would essentially only lead to an additional drift of the qubit state, i.e., a ran- 
dom walk with a nonzero average value. Neither our analytical results nor our conclusions 
would otherwise change qualitatively. In fact, assuming the the bang-bang pulse cycles are 
sufficiently short relative both the typical +1 to —1 and —1 to +1 switching times of the 
bfl, bang-bang suppression of the bfl noise should not be diminished at all by the bfl's asym- 
metrical switching. We therefore obtain Equ. as our starting point of the bfl-perturbed 
qubit dynamics. 
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